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1 Introduction 

Quantum chromodynamics (QCD) dH[2]) is the accepted theory of strong nuclear 
forces. But QCD is hard to solve. Perturbation theory gives access to physics 
at high scales, well above the confinement transition. Lattice gauge theory is 
best suited for the calculation of static quantities such as the vacuum structure, 
the spectrum, and thermodynamics. Various heuristic models and approximation 
schemes, ranging from the bag model to fluid dynamics, supply partial informa- 
tion about the dynamics of QCD in the strongly coupled regime. 

Over the past ten years or so, string theory has supplied a new way of thinking 
about strongly coupled gauge theories: the gauge-string duality ([31H1E]). One 
of the simplest cases is maximally supersymmetric Yang-Mills theory, which is 
understood to be dual to string theory in a specific curved ten-dimensional ge- 
ometry, AdS§ x 5 5 . The gauge-string duality can also handle confining gauge 
theories which are at least qualitatively similar to QCD. 

As compared to the more standard tools for understanding QCD, the gauge- 
string duality offers two main advantages. First, it offers direct access to the 
strong coupling region. The simplest examples, like maximally supersymmetric 
Yang-Mills theory, are perfectly conformal, with identically zero beta function. 
But they are nevertheless strongly interacting, and confining behavior is in some 
sense just around the corner. Second, it is straightforward to introduce finite 
temperature without passing to Euclidean signature: one instead introduces a 
black hole horizon. This makes it relatively straightforward to extract transport 
coefficients, such as viscosity and diffusion constants, for analogs of the quark- 
gluon plasma. All these real time phenomena are challenging to calculate on the 
lattice; the complex weighting factor e lS where S is the action in the path integral 
makes the standard importance sampling used on the lattice unreliable. 

The gauge-string duality also suffers from several deficiencies as a phenomeno- 
logical tool. The worst one is that we do not have a dual for QCD itself. We are 
therefore forced to make comparisons in a more or less ad hoc fashion between 
real-world QCD and the theories we can deal directly with via the gauge-string 
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duality. Another serious problem is that the strong coupling limit is taken quite 
strictly in many or most gauge-string calculations: only the leading behavior in 
the limit of infinite coupling is readily accessible. Significant inroads have been 
made on computing subleading corrections, but both technical and interpreta- 
tional challenges are substantial. Likewise, it is in general difficult to proceed 
past the leading order in large N expansion, where iV is the number of colors. 

The role of the gauge-string duality in understanding strong interactions is 
still evolving. The duality has made an indelible mark on our understanding 
of QCD, and it offers uniquely geometrical perspectives on confinement and the 
quark-gluon plasma. But it hasn't solved QCD. 

The rest of this review is organized as follows. In section [2] we introduce the 
basic ideas of the gauge-string duality. In section [3] we summarize how flavor 
is included by introducing branes into negatively curved bulk geometries in the 
probe approximation. In section [5] we outline gauge-string constuctions that 
capture aspects of confinement. In section [5] we present gauge-string calculations 
of shear and bulk viscosity of strongly coupled, deconfined plasmas. In section [6] 
we review work on expanding plasmas in the dual supergravity description. In 
section [7] we summarize several approaches to the physics of hard probes of a 
thermal medium. We conclude briefly in section [8l 

The size of the relevant literature and the brevity of this review have prevented 
us from discussing many of the topics with the thoroughness they deserve. Also, 
we have altogether omitted some important developments, notably the work on 
anomalous dimensions of long operators (see for example ([6])); string descriptions 
of gluon scattering ((7J); bottom- up approaches to holographic models of confine- 
ment (see for example (8,9,10)); and light-front holography, reviewed for example 
in CLU). 

2 The basics of the gauge-string duality 

The story of the gauge-string duality turns on understanding D3-branes. There 
are other types of branes in string theory, but D-branes are the best understood, 
thanks to the crucial contribution of Polchinski (|12j) . By definition, D-branes are 
locations in ten-dimensional spacetime where strings can end. A D3-brane is a 
D-brane extending in three spatial dimensions as well as in the time direction. 
The low-energy excitations of a single D3-brane are described by N = 4 super- 
symmetric £7(1) gauge theory. If N D3-branes are placed on top of one another, 
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one finds instead N = 4 supersymmetric U(N) Yang-Mills gauge theory (|13p 
(hereafter M = 4 SYM). This theory splits into a U{\) part, which is free, and 
relates to the center of mass motions of the stack of D3-branes; and an SU (N) 
part, which is interacting, and relates to the relative motions of the branes. The 
gluons of these gauge theories can be represented as strings running from one 
brane in the stack to another. If N = 3 and we label the branes as R, G, and B, 
then a string running from R to B has color quantum numbers RB, as appropri- 
ate for a gluon of SU (3) gauge theory. The M = 4 supersymmetry of D3-branes 
dictates that this RB string can act not only as a gluon, but also as spin-1/2 
fermion or as a scalar; the strings describing these superpartners of the gluons 
differ from the ones describing the gluons themselves by fermionic excitations on 
the string. Regardless of spin, all particles in N = 4 SYM are charged in the 
adjoint of the gauge group. The precise dynamics of the gauge theory is almost 
entirely dictated, at the level of a renormalizable lagrangian, by the number of 
branes N and the supersymmetry. The only freedoms one has are to adjust the 
gauge coupling and a 9 parameter. The 9 parameter breaks CP, and we set it to 
zero from here on. This amounts in the D3-brane construction to requiring that 
a particular scalar field (more precisely, a pseudoscalar) should vanish. 

D3-branes have a definite mass per unit volume, and they have a definite 
charge under a 5- form field strength. They therefore deform spacetime into a 
solution of the ten-dimensional Einstein equations coupled to the 5-form. Close 
to the D3-branes, this solution takes the form of a direct product, AdS^, x S 15 , 
where AdS stands for anti-de Sitter space. The gauge-string duality ([UIUH]) 1S 
an equivalence between string theory on this geometry and J\f = 4 SYM in four 
dimensions. Green's functions of gauge-invariant operators, Wilson loops (|14pi5p . 
and thermal states (|16[I17|) are just some of the gauge theory constructions which 
have well-studied translations into the language of string theory in anti-de Sitter 
space. 

To get started with AdS/CFT, the first step is to understand the metric of 
extremal D3-branes. Extremal in this context means that the D3-branes are at 
zero temperature, so they are described by a vacuum state of Af = 4 super- Yang- 
Mills theory. The metric is 



ds 2 w = H-Wi-dt 2 + dx 2 ) + R x l 2 (dr 2 + r 2 dVL 2 ) , 



(1) 



where 




(2) 
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and L is related to the number of D3-branes: 

LA = ^ = &mN^ • (3) 

Here n = \/8ttGiq is the gravitational coupling in ten-dimensional supergravity, 
9ym = \J 4vrg str ing is the gauge coupling of N = 4 super- Yang- Mills, N as before 
is the number of D3-branes, and a 1 is the Regge slope parameter of fundamental 
strings, inversely proportional to their tension. Often we will use the notation 
A = gyjtfN for the 't Hooft coupling. 

The region where r S> L is asymptotically flat ten-dimensional spacetime, R 9,1 . 
The D3-branes can be thought of as a defect in that spacetime, extended in the 
t and x directions, and having a characteristic size L in the other six directions. 

The crucial technical step of ([3]) is to focus on the region r C L by dropping 
the 1 from the function H. The resulting geometry is a direct product of AdS§ 
and S 5 : ^ 2 

ds\ Q = ^{-dt 2 + dx 2 ) + ^jdr 2 + L 2 dQ 2 . (4) 

Li T 

Physically, dropping the 1 from H means we're zooming in close to the D3-branes: 
so close that we lose track of the asymptotically flat region. In short, is a 
way of replacing the branes by a curved geometry, and the essential claim of 
AdS/CFT is that the gauge theory dynamics built from strings on the branes is 
equivalently captured by the geometry (J4j). 

A variety of different choices of radial variable are in common use. Among the 
most common are the following two: 

L 2 r 

Z=— U =j2- ( 5 ) 

Note that z has dimensions of length, while u has dimensions of energy. Using 
the z coordinate makes it clear that AdS§ is conformal to fiat space: 

L 2 

dsl = —(-dt 2 + dx 2 + dz 2 ) . (6) 

More precisely, AdS§ is conformal to the z > half of five-dimensional flat space, 
R 4,1 . It has a boundary at z = 0, and that boundary is Minkowski space, R 3 ' 1 , 
parametrized by t and x. A common way of speaking is that the gauge theory 
lives on this boundary, while string theory lives in the bulk of AdS§. This point 
of view is useful, but it is perhaps more accurate to say that string theory on the 
whole of AdS& x <5 5 is equivalent to M = 4 super- Yang-Mills on R 3 ' 1 . 

To get a feel for the geometry of AdS§, consider the following two thought 
experiments: 
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1. An observer at fixed depth zq shoots a light ray straight upward toward 
the boundary. What happens? Well, the light ray travels upward along a 
trajectory z(t) = Zq — t, so it reaches the boundary in a time t = z$. At 
least some of the light presumably is reflected back down to the observer, 
who receives it at a time t = 2zq. This is all exactly as it would be if the 
observer lived in the z > half of R 4,1 . 

2. The observer extends a rod upward toward the boundary of AdS§. His 
rod is very long. He notices however that its tip never hits the boundary. 
What's going on? Well, the proper distance between the observer and the 
boundary is infinite, because it is given by the integral 



The divergence comes from the region near z = 0. Evidently, this is very 
unlike the situation in the z > half of R ' . 

In summary, AdS§ gets very, very big near the boundary, but light can neverthe- 
less get to the boundary in finite time. 

A crucial ingredient in the relation of AdS§ x S 5 to J\f = 4 SYM is the action of 
dilations. First, on the field theory side, the beta function of M = 4 SYM vanishes 
identically. So the renormalization group says that dilations act trivially on the 
theory: T m m = as an operator identity. If we turn to the gravity side, a dilation 
should be just a coordinate transformation, x m — > Kx m , where x m = (t,x,y,z) 
and K is a constant. By itself, this coordinate transformation does not preserve 
the form of the metric Q. But if we also send r — > r/K, or equivalently u — > u/K 
or z — > Kz, then the metric is preserved. Thus we have associated dilation 
symmetry of N = 4 SYM with an isometry of the AdS§ metric. Taking K > 1 
means making things bigger in the x m directions. At the same time, u — ► u/K 
means that the location in the radial direction gets smaller — "deeper" in AdS^, 
and further from the boundary. 

The upshot of this discussion is that the region of AdS^ where u is large cor- 
responds to ultraviolet (UV) physics, whereas the region where u is small corre- 
sponds to infrared (IR) physics. A slight generalization of the AdS part of the 
metric (jlj) helps clarify this: 




(7) 



ds 2 = L 2 u 2 (-hdt 2 + dx\ + dx\ + dxj) + — + L 2 dttj 



(8) 



where h = 1 — u^/u 4 . This geometry satisfies the Einstein equations with the 
same five-form field strength that supported the original AdS§ x 5 5 metric. The 
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solution (jHJ) describes the near-horizon geometry of near-extremal D3-branes. 
Near-extremal means that a little bit of mass has been added to the D3-branes, 
in the form of thermal energy. The temperature manifests itself as the Hawking 
temperature of the horizon at u = Uh, which is T = u^/vr. Thus we can regard u 
or u/tt as a typical energy scale pertaining to processes taking place in AdS§ at 
a depth u. 

To understand the calculations of Green's functions in the gauge-string duality, 
it helps to return for a moment to the idea of a D3-brane as a defect with 
transverse size L in asymptotically flat ten-dimensional spacetime. In order to 
explore its physical properties, the first thing we might try is throwing things 
at it and seeing how they bounce off or are absorbed. If not for the fact that 
the D3-brane extends infinitely far in the three directions parametrized by x, 
we could assign cross-sections to processes of this type. Because of its infinite 
extent, the right quantity to consider is the cross-section per unit volume in the x 
directions. There is a large literature devoted to exploring D-branes in this way, 
including (18 HT9l[20| l2i p22U23ll . and the gauge-string correspondence is in part an 
outgrowth of it. The typical process by which a D3-brane absorbs a closed string 
is to split it into two or more open strings that propagate along the worldvolume 
of the branes. If the closed string is a graviton, such a process is mediated by 
the coupling h mn T mn , where T mn is the stress tensor of M = 4 SYM. Then the 
open strings are gluons or superpartners of gluons. From the field theory point 
of view, the calculation of the response of the brane to the perturbation by an 
external metric amounts to the calculation of correlation functions of the stress 
tensor. 

To compute the absorption cross-section of a graviton by D3-branes, one is 
interested in an inclusive process where the final state is arbitrary. The optical 
theorem guarantees a relation between the cross-section of the inclusive process 
and the two-point function of the vertex operator T mn . The two-point func- 
tion characterizes the linear response of the D3-branes to a metric perturbation. 
Higher point functions would characterize this response past the linearized ap- 
proximation. The calculation of Green's functions dUE]) can be understood, in 
part, as a streamlined treatment of absorption processes. Here is a heuristic 
summary how gauge theory Green's functions are calculated using ^4<i55 x5 5 . A 
graviton perturbation h mn creates an environment that the D3-branes respond 
to. But in the gravitational picture, the near- horizon geometry (i.e. AdS§ x S 5 ) 
effectively replaces the D3-branes. So the graviton perturbation creates an envi- 
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ronment that AdS§ x S 5 responds to. That environment should be thought of as 
imposed at some very large value of u, because that is where AdS§ x S 5 matches 
onto asymptotically flat space. So the prescription should be to impose some 
metric deformation in the limit of large u and then carry out the path integral of 
string theory in ^4<i55 x S 5 subject to boundary conditions corresponding to the 
deformation. The result of the path integration then has to be the generating 
functional of Green's functions for color-singlet operators in the field theory. If 
the metric is the only field one deforms at large u, then one has access only to 
n-point functions of T mn . But there is a large list of additional operators for 
which dual fields in AdS§ x S 5 are well understood. 

A big stumbling block is that the path integral of string theory in AdS§ x 5 5 
is poorly understood. The best-studied approximation to it is the saddle point 
expression e lS , where S is the gravitational action, classically extremized subject 
to appropriate boundary conditions. This approximation is good when AdS§ has 
small curvature compared both to the Planck scale and to the string scale. The 
former constraint turns out to mean that the number of D3-branes must be large: 
N 3> 1. The latter means that the 't Hooft coupling must be large: gyM-^ ^ 1- 

3 The addition of fundamental matter 

One ingredient that our discussion so far has been missing is fundamental matter. 
The quarks of QCD are fermions transforming in the fundamental representation 
of the SU(3) gauge group. They interact via exchange of gluons, which are 
massless spin-1 bosons transforming in the adjoint representation of SU(3). The 
J\f = 4 theory describes supersymmetric glue, in which all particles transform 
in the adjoint. So we are still missing the quarks. In the limit of a large num- 
ber of colors that underlies the gauge/string correspondence, the difference in 
representation is crucial. Topologically, the fundamental representation quarks 
introduce boundaries into the double-line Feynman diagrams. Each fundamental 
matter loop suppresses the diagram by a power of N in the large N counting. 
For example, the free energy of a free gas of gluons and quarks has an order iV 2 
contribution from the gluons and an order NfN contribution from the quarks, 
where Nj is the number of different flavors of quarks. The naive large N counting 
will work as long as Nf <C N. For this very basic N counting the spin of the 
quarks is irrelevant; for that reason it has become common in the gauge/string 
literature to refer to all fields in the adjoint representation collectively as "the 
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glue" and all fundamental representation fields as "the quarks" or "the flavors", 
irrespective of spin. One important consequence of this large N counting is that 
dynamical fermion loops are suppressed at large N. Inclusion of the flavor fields 
allows for new observables, e.g. the chiral condensate, the meson spectrum, flavor 
current correlators or the order iV contribution to the free energy. But all inter- 
nal lines in the diagrams are only glue lines. The quark dynamics is happening 
in the background of the strongly coupled glue. One can for example even add 
additional fundamental matter to the conformal J\f = 4 SYM. This will make 
the [3 function positive and asymptotic freedom is lost. The coefficient of the (5 
function however is of order 1/N and so this conformal symmetry breaking and 
the appearance of the Landau pole can be ignored at large N. 

These statements about Feynman graphs have of course also a direct analog 
on the string theory side. Since in the topological counting fundamental matter 
introduces boundaries, the string theory dual will have to have open strings, 
too. Including open strings amounts to the addition of D-branes in the bulk, the 
places where strings can end. D-branes in string theory have a tension of order 
one over the coupling, that is of order N in the limit of large TV at finite g s N; 
the stress tensor xj^ hra,ne hence also is of order N. Newton's constant Gn on 
the other hand is of order g§ ~ A -2 . In the presence of D-branes Einstein's 
equations need to be re-solved as the additional D-branes curve space. However 
the right hand side is of order Gn Tj/*~ hra " nc ~ 1/N and hence vanishes in the 
large N limit. The classical dynamics of these "flavor D-branes" will demand 
that they minimize their action in the presence of a given background geometry. 
This essentially boils down to a minimal area problem, as the action of a D-brane 
is proportional to its worldvolume. (In some situations it is important to know 
that there is a gauge field on the brane, unrelated to color. Its inclusion leads 
to the Dirac-Born-Infeld action or DBI action. There are also fermionic fields on 
the brane.) As in the field theory the dynamics of the glue, which is encoded 
in the background geometry, is unaltered by the presence of the flavor branes at 
leading order in N. The first example of such a flavored version of AdS/CFT 
was a conformal example with massless quarks but with the backreaction fully 
included (|24p . The realization (|25U26j) that flavors, massless and massive, can 
be studied in the "probe limit," that is in the limit where the backreaction of 
the brane can be neglected, opened the door to a systematic study of the flavor 
sector in AdS/CFT. 

In order to make these ideas more explicit, let's consider in more detail a D7- 
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brane probe in AdS§ x S 5 , following (|26p . It is useful to begin with what seems 
like a detour: consider the D7-brane and the D3-branes both in the probe approx- 
imation in ten-dimensional flat space. Let's assume that the D7-brane is extended 
along the same spatial directions as the D3-branes, and that the distance of clos- 
est approach between the D7 and the D3's is ro7- It is a textbook exercise to 
show that eight supersymmetries are preserved by such a configuration — provided 
we leave all the branes at zero temperature. The dynamics of low-energy exci- 
tations, including strings running between the D3-branes and the D7-brane, is 
completely determined by the supersymmetry. The D3-D7 strings, in particular, 
are described in terms of massive scalars and fermions, fundamentally charged 
under the SU(N) gauge group on the D3-branes. Their mass is m = r^flna 1 , 
simply because this is the tension of a string with length td7- 

With fully a quarter of supersymmetry preserved, one expects that the behavior 
of the D7-branes will scarcely be altered if one replaces the probe D3-branes by 
the warped metric ([I]). The dr 2 + r 2 dQ, 2 part of this metric is the six-dimensional 
flat space transverse to the D3-branes. The D3-branes are located at the origin 
of this transverse space. The simplest embedding of the D7-branes is to have 
them run along a four-dimensional plane in the transverse space whose closest 
approach to the origin is at a distance r^j. If we express the metric of the unit 
S 5 as 

dO§ = dO 2 + cos 2 dip 2 + sin 2 dUj , (9) 

where dQ, 2 is the metric of a unit S* 3 , then the four-dimensional plane can be 
specified by the equations ip = 0, rcosO = td7- It can be checked that this em- 
bedding does extremize the D7-brane worldvolume action, which is proportional 
to its area. The mass of a string running from the D3-branes to the D7-brane is 
still m = rD7/27ra', because in the string world- volume action, a factor of i? 1//4 
from the string's extension in a transverse direction cancels against a factor of 
jj— 1/4 f rom j£ g propagation in the time direction. 

To achieve a description of a D7-brane probe in AdS§ x S 5 , one can rely upon 
the previous discussion, but simply drop the 1 from the harmonic function H. 
Geometrically, the D7-brane wraps an S 3 in the <S 5 , and it extends across a 
range of radii in AdS§, r > rj^-j. A string stretched between the D3-branes to 
the D7-brane is now understood as a string running from the depths of AdS§, 
where r — > 0, up to r = r^j. One may also consider as a limiting case r-^j = 0, 
which means that the D7-brane passes right through the D3-branes. Then the 
strings running from the D3-branes to the D7-brane act as massless fundamentally 
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charged particles. In this special case, the D7-brane wraps an equatorial S 3 in 
5 5 . It might appear that an instability should arise from perturbations where the 
D7-brane slips away from the equator to an S 3 of smaller size. In fact, this mode 
corresponds to a tachyon of mass m = -3/L 2 , but according to the analysis of 
(127h . tachyons do not indicate an instability in AdS§ until m 2 falls below —4/L 2 . 
Heuristically, this is because AdS§ acts as a box, and fluctuations must always go 
to zero at its boundary. The lowering of potential energy one achieves by exciting 
a tachyonic mode is offset, provided m 2 > —4/L 2 , by the raising of kinetic energy 
due to this boundary condition. 

The complete meson spectrum in the example of the D3/D7 system was worked 
out in (I28p . Mesons are the eigenmodes of small gauge field fluctuations of the 
brane. Gauge field propagation on the brane in response to an excitation by an 
external global current calculates current correlation functions in the dual field 
theory. This is in direct analogy with the relation between n-point functions 
of the stress tensor in the field theory and the analysis of the response of the 
AdS$ geometry to an external gravitational perturbation we described before. 
For massive quarks this correlation function is dominated by discrete modes, the 
mesons. 

Other types of flavor can be constructed by starting with D-branes of different 
dimensions. Let there be N coincident Dp-branes, which we refer to as "color 
branes" because the gauge theory on their world-volume is thought of as the 
one analogous to SU (3) color. Let there be one (or Nf <C N) D^-branes, called 
"flavor branes," which are treated in the probe limit. The field theory description 
of such a setup depends crucially on one parameter: the number of so-called "ND" 
directions. An ND direction is a direction which is parallel to the world- volume 
of the flavor branes, but orthogonal to the world- volume of the color branes, or 
vice versa. In the D3/D7 example above, there are four ND directions, because 
the D7-branes extend in four directions that the D3-branes do not. Systems with 
four ND directions preserve eight supersymmetries, which essentially fixes the 
Lagrangian. Another interesting class of flavor branes has six ND directions. In 
this case all supersymmetry is broken. The fundamental matter that is added is 
just a single Weyl fermion. We will discuss one example of this type below, the 
Sakai-Sugimoto model. 
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4 Gravity Duals of Confining Gauge Theories 



4.1 The Hard Wall paradigm 

As scale transformations map to an isometry of AdS$, it is clear that in order 
to describe a confining gauge theory we need to move beyond AdS§ in the dual 
gravity. More particularly, energy scale in the field theory is encoded in the radial 
direction of the dual gravitational theory, so there will be a spatial radial coordi- 
nate z* that maps to the strong coupling scale Aqcd- in pure AdS§, plane- wave 
normalizable fluctuations around the background can come with any timelike 
four- momentum p 2 . These modes should be thought of as the particle spectrum 
of the dual field theory. This continuous spectrum is completely appropriate for 
the dual of a conformal field theory without well defined particles. For a confining 
gauge theory with a discrete excitation spectrum, the dual gravity side should 
also only have discrete modes for the fluctuations with four momentum p 2 = 
with discrete masses m n . This can easily be achieved if the extra dimension 
truncates around instead of extending all the way to z = oo. Only a discrete 
set of modes will satisfy both the boundary condition at z = as well as at z*. 
The details of the confinement mechanism are encoded in the geometry around 
z*. Many generic features of confining theories with a gravity dual do not depend 
on these details and can be analyzed by simply truncating undeformed AdS$ at 
z* as first advocated in (|29p . This is often referred to as a "hard wall." In the 
following, we will review several string theory constructions that realize this hard 
wall paradigm. In most of the known examples the geometry is actually quite 
distinct from AdS§ even in the UV. Only one of the examples we will discuss 
becomes AdS§ in the UV. 

Despite the fact that these are confining gauge theories, the confining theories 
with gravity dual are still very different from QCD. In the original conformal 
example of AdS / CFT we need to work in the limit of large 't Hooft coupling in 
order to be allowed to approximate the classical bulk string theory by classical 
supergravity. One way to understand this is that large 't Hooft coupling ensures 
a large hierarchy between the mass of the excited modes of the string (governed 
by the string length) and the light modes of the string (governed by the curvature 
radius) . The confining theories also have to have this large hierarchy of scales in 
order for classical gravity calculations to be applicable. Again the string scale on 
the gravity side (dual to the flux tube tension in the field theory) is much larger 
than the typical energies of the gravitational modes (which map to the mass gap 
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of the field theory). A typical confining gauge theory will have a flux tube tension 
a set by the same scale Aqcd as the mass gap. In short, a confining theory with 
a gravity dual needs a large dimensionless number A ~ a 2 /Aq CD . 

Confining theories with gravity dual share many common qualitative features 
which can be understood in terms of the hard wall paradigm. Most importantly, 
Wilson loops will exhibit area law behavior. To calculate the Wilson loop one 
needs to identify the minimal area worldsheet connecting the two test charges at 
the AdS boundary (|15yi4|) . As for all other correlation functions, this worldsheet 
is the response of the gravity system to a perturbation imposed from the out- 
side: in this case a quark and anti-quark pair that remain at fixed separation for 
a specified time. In the conformal case, the worldsheet sinks down deeper and 
deeper into the bulk as the external quarks are separated along the field theory 
directions. In a theory with any kind of hard wall, the string can sink no further 
than z*. For large separation between the test quarks the minimal area is hence 
dominated by a string lying at the bottom of the geometry at z* . The latter cor- 
responds to a flux tube with an effective tension given by \J gu (-z* ) 9xx ( z * ) /27ra' , 
where l/27ra' is the tension of a fundamental string. For N = 4 SYM with a 
hard wall at = 1/A this evaluates to ^A 2 . Another property confining gauge 
theories with gravity dual share with their generic cousins is power laws in elastic 
hard scattering processes (29) and Regge physics (|30|) . that is scattering at large 
center of mass energy s at fixed Mandelstam parameter t. 

One place were the difference between generic confining gauge theory and a 
confining gauge theory with gravity dual becomes apparent is deep inelastic scat- 
tering (DIS) as studied in (|3Tj) . DIS requires one to understand the time-ordered 
two point function of an electromagnetic current in an hadronic state. At large 
q 2 the matrix element can effectively be calculated using the operator product 
expansion (OPE). So the main physical question is what operators dominate the 
current-current OPE. At weak coupling it is simple single trace operators that 
dominate the OPE. In QCD they are bilinear in the quark fields. They create 
and then annihilate a single parton from the vacuum. At strong coupling these 
operators generically get an anomalous dimension of order A 1 / 4 . In the gravita- 
tional dual the fields corresponding to these operators are excited string modes. 
The dominant operators must be dual to supergravity modes, not string modes. 
In (|3ip it was argued that the corresponding field theory operators are double 
trace operators, that is they are bilinear in hadron fields, and hence create and 
annihilate a whole hadron, not just a single constituent. A more intuitive way 



14 Gubser and Karch 

to describe this result is in terms of the parton content of a hadron. At strong 
coupling, that is for confining gauge theories with a gravity dual, processes are 
typically dominated by a very large number of very soft gluons. This stands in 
stark contrast to a weakly coupled theory such as QCD in the asymptotically 
free regime where, in a typical process, a few partons carry order 1 fractions of 
the total energy and momentum. This basic qualitative insight guides the under- 
standing of many of the strong coupling processes at zero and finite temperature 
we will describe in this review. Since soft gluons are the main source of diffi- 
culty in QCD, gauge theories with supergravity duals that are dominated by soft 
gluons are a great theoretical laboratory. 

4.2 Mass Deformed M = 4 theories: J\f — 2* and N = 1* 

Closest among confining backgrounds to the hard wall paradigm described above, 
with a spacetime that is AdS§ in the UV and ends in a wall in the IR, is the 
Polchinski-Strassler (PS) background (32). On the field theory side the idea is 
extremely simple: in order to reduce N = 4 SYM to pure glue, add mass terms 
to the extra scalars and fermions. In the UV the mass terms are inconsequential, 
so asymptotically the dual spacetime is the standard AdS^ x S* 5 solution of type 
IIB string theory. In the IR the heavy fields can be integrated out and one is left 
with pure glue, which is expected to be a confining theory. Flavor can be added 
using the technique described in section [31 as was done for the PS background in 

AMD- 

How about the hierarchy between flux tube tension and mass gap, which we 
argued above should be a feature of any confining theory with a gravity dual? 
When we say that at scales below the mass of the fermions and scalars the theory 
reduces to pure glue, we have implicitly been assuming that we can take a limit 
where the mass m is much larger than Aqcd- So what happens if we take this 
limit m 3> Aqcd? At scales below m the theory is pure glue and we know how 
the coupling evolves with scale. Asymptotic freedom of the pure glue theory 
guarantees that by the time we let the coupling evolve from Aqcd tom> Aqcd 
the Yang-Mills coupling is small. At energies above m we have the conformal 
J\f = 4 theory and the coupling does not evolve any further. So the limit of 
m 3> Aqcd forces us to consider the M = 4 theory at weak coupling. But this is 
not where the supergravity analysis is reliable! Conversely, if one sticks with the 
strongly coupled N = 4 at energies above m so that supergravity is reliable, one 
is forced to take m ~ Aqcd- The theory one obtains this way does not reduce 
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to pure glue: it has extra light adjoint fields with masses of order Aqcd- These 
extra light fields conspire to push the flux tube tension up to a parametrically 
higher scale than the scale of the mass gap. It is, of course, also possible to 
only give mass to some of the extra fields leaving one with pure super Yang- Mills 
in the limit m S> Aqcd 5 but again with extra light adjoints in the limit that 
supergravity is valid. The corresponding gauge theories are typically referred to 
as J\f = 2* and N = 1*, depending on the amount of unbroken super symmetry. 
The superscript * is understood to indicate the presence of the extra light adjoint 
fields. 

4.3 Witten's Black Hole and Sakai-Sugimoto Model 

One simple way to break both supersymmetry and conformal invariance is to 
heat a system up to finite temperature. In Euclidean signature, this corresponds 
to making time periodic with anti-periodic boundary conditions for the fermions. 
Long distance questions are dominated by an effective theory, which for both 
QCD and J\f = 4 SYM is just pure Yang-Mills theory in 3 real dimensions: the 
fermions are massive due to the anti-periodic boundary conditions, the scalars 
inherit this mass via loop corrections and the electric field is heavy due to De- 
bye screening which gives the At component a mass of order gT. Long dis- 
tance finite temperature physics is the physics of magnetic fluctuations. The 
fact that pure glue in 3d confines makes long distance physics at scales l/(g 2 T) 
non-perturbative even at arbitrary high temperatures in QCD. One can reinter- 
pret these well known facts to achieve the goal of realizing a confining gauge 
theory with gravitational dual: instead of viewing this system as N = 4 SYM 
or QCD at finite temperature one can also view it as a 3d gauge theory at zero 
temperature which was obtained by compactifying N = 4 SYM on a spatial cir- 
cle of radius R = ^ with anti-periodic boundary conditions for the fermions. 
At low energies this 3d theory is a theory of pure glue with a dynamical scale 
A ~ 9ym 3d = 9ym 4d(R)/R- There are additional degrees of freedom which come 
from the higher Fourier modes of the 4d fields on the circle, the so-called Kaluza- 
Klein (KK) modes. They have a mass of order 1/R. For a weakly coupled theory 
(which in the case of QCD means very small R) the dynamical scale and the mass 
scale of the KK- modes are widely separated. In the limit that the 4d coupling 
and R both go to zero keeping the 3d coupling fixed, the KK-modes completely 
decouple, and in this limit the theory describes pure 3d glue. It is the opposite 
limit of strong coupling in 4d in which the theory has a gravity dual. In this limit 
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we have a theory of 3d glue coupled to the tower of KK-modes whose mass is of 
the same order as the dynamical scale. Again, these additional modes conspire 
to push the flux tube tension up to a much higher scale than the mass gap. 

As we already discussed, in the gravity dual finite temperature corresponds 
to the introduction of a black hole; the metric is given by the Wick rotation of 
Eq. ([8]). The Euclidean black hole geometry relevant for equilibrium physics 
has the geometry of a cigar: Euclidean time is a circle, and it shrinks to zero 
size at the horizon, the tip of the cigar. One may reinterpret this geometry as 
a 2+1 dimensional confining gauge theory by Wick rotating one of the spatial 
coordinates back to be the time direction of this 2+1 dimensional gauge theory. 
Then the metric in Eq. ([8]) becomes 

u 2 (-hdt 2 + dx 2 + dy 2 + dz 2 ) + ^ - u 2 (hdr 2 + dx 2 + dy 2 - dt 2 ) + (10) 

The same strategy can also be used to construct a gravity dual for a 3+1 dimen- 
sional confining gauge theory (|17p . One simply reinterprets the gravity dual of 
maximally supersymmetric YM theory in 4+1 dimensions at finite temperature, 
which is given in terms of the near horizon geometry of near-extremal D4 branes 
in type IIA string theory, as a confining gauge theory in 3+1 dimensions. The 
gauge theory has a tower of KK-modes in addition to a theory of pure glue in 
3+1 dimensions. In the weak coupling limit one describes pure glue; in the limit 
where the gravity dual is weakly curved, the KK-modes have masses of order 
Aqcd- The resulting metric is (i, j run over the 3 spatial directions, xt± is the 
compact circle to get from 5 to 4 dimensions) 



ds 2 



Rda) 
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3/2 



du 

— — + u^d\l A 



-dt 2 + bijdx % dx 3 + f(u)dxj 

(11) 

where f(u) = l-(^) and ua is the minimum value of the u-coordinate, the tip 
of the cigar. There is also a scalar field turned on: the 5d dilaton has a profile 

, / \ 3/4 
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In the field theory, R is the radius of the circle one compactified 4+1 SYM on, 
g 2 = (2ir) 2 g s l s is the 5d YM coupling, and correspondingly g\ = g 2 /(2irR) is the 
4d YM coupling. 

To add fundamental flavors to this geometry, we can follow the approach de- 
scribed in Section [3l Two obvious candidates are probe D6-branes (f34"|) or probe 
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D8-branes (|35j) with four and six ND directions respectively. The former pre- 
serves eight supersymmetries of the 5d super Yang-Mills. But supersymmetry in 
the flavor sector is broken by the compactification with anti-periodic boundary 
conditions for the fermions, just as it was in the glue sector. Chiral symmetry 
ensures massless fermions in 4d, but the scalar superpartners become massive via 
loop corrections. In the weak coupling limit one obtains real QCD with fermionic 
quarks only. In the strong coupling limit one has even more extra matter with 
masses of order the strong coupling scale: the KK-modes of the flavor fermions 
and their scalar partners. The advantage of using the D8-brane probe, with six 
ND directions, is that one introduces only fermions. A single D8 introduces a 
chiral fermion. Anomaly freedom requires that the D8 is accompanied by an 
anti-D8. This system is stable as long as the D8's are separated in the direction 
of the compactified circle. 

Both of the systems described in f)34|.l35j) nicely exhibit chiral symmetry break- 
ing of the axial U(l) symmetry (whose anomaly is negligible at large N). The 
low energy excitations of the flavor branes are the mesons (I28p . and hence the 
DBI action plays the role of the chiral Lagrangian. The D8 model has the dis- 
advantage that it is difficult to add a mass term for the fermions. As fermions 
of different chirality live at different places in the fifth dimension, the mass op- 
erator is non-local. To some extent the effects of such a non-local operator can 
be captured via the standard dictionary (|36p . The advantage of the D8 system 
is that for more than one flavor one can observe chiral symmetry breaking of 
the full non-abelian chiral SU(Nf) x SU(Nf) to its diagonal subgroup in a very 
beautiful geometric manner: deep in the IR of the geometry the D8 and anti-D8 
reconnect. Due to confinement, the D8's cannot end in the IR and have no op- 
tion other than to connect. Thus one sees in this model that confinement implies 
chiral symmetry breaking. The D8 system also has an additional parameter, 
the separation of the 4d flavors along the compactification circle. This allows 
one to dial the physics in the meson sector independently of the glue sector. In 
particular, the typical masses in the meson sector are set by L, the asymptotic 
separation of the D8 and anti-D8, while the glueball masses are set by R, the 
compactification radius. For 1/L 3> 1/R the mesons are much heavier than the 
glueballs. For 1/L S> 1/R one also has a very interesting phase structure where 
in the temperature window between 1/R and 1/L the theory is deconfined, but 
chiral symmetry is still broken. The D8 system, often referred to as holographic 
QCD (HQCD) or the Sakai-Sugimoto model, is currently the closest to QCD we 
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have come with a theory that has a gravity dual. 
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4.4 The warped deformed conifold or Klebanov-Strassler Back- 
ground 

Another popular supergravity solution dual to a confining gauge theory is the 
warped deformed conifold of Klebanov and Strassler ([37|) . The dual gauge theory 
is super symmetric and exhibits (discrete) chiral symmetry breaking even before 
adding additional flavors via probe branes. The UV of this theory is very inter- 
esting: the dual field theory has a product gauge group out of which one factor 
is asymptotically free, while the other runs to strong coupling in the UV. As this 
second gauge factor approaches its Landau pole, the gauge bosons turn out to be 
better described as composites of a different gauge group. Using these dual vari- 
ables the roles of asymptotically free and non-free gauge group factors change and 
one can run further up in energy until the other gauge coupling becomes strong 
and the same story repeats itself. Like in our other examples, one can take a 
limit in which this theory reduces to pure glue in the IR, but this involves string 
scale curvatures in the bulk. The limit in which supergravity is applicable has 
again extra matter of mass of order Aqcd whose net effect is to push up the flux 
tube tension. Reviews have been written about the warped deformed conifold 
and its applications to field theory and cosmology: see for example (|38|,l39p . 

4.5 Confinement/Deconfinement transition 

One more universal feature of confining theories with a gravity dual which we 
haven't yet discussed in detail is their phase structure as a function of tempera- 
ture. Three-color QCD with 2 light flavors and a strange quark can have a first 
order transition, a second order transition, or a smooth crossover from hadronic 
degrees of freedom at low temperature to a gas of quarks and gluons at high tem- 
peratures depending on the precise values of the quark masses (with real world 
QCD exhibiting a crossover according to lattice simulations). In contrast to this, 
confining large A glue theories are expected to always have a first order transition 
separating the hadronic phase from the quark-gluon phase. The reason for this 
is that the A scaling of the free energy is very different in the two phases. The 
confined hadrons have order 1 degrees of freedom, whereas the deconfined glue 
has order A 2 . The order A 2 part of the free energy is identically zero at low 
temperatures, so that the derivative of the free energy has to be discontinuous at 
the critical temperature T c where we go over from hadron gas to gluon plasma. 
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The presence of fundamental matter does not alter this conclusion as long as the 
number of flavors Nf is finite. In the deconfined phase, the quarks contribute a 
sub leading order NfN contribution to the free energy. Only when Nf ~ N can 
one get a smooth crossover as in the real world. In this case the Nf mesons in 
the confined phase can smoothly connect with the N 2 gluons in the deconfined 
phase. 

In the gravity dual, this phase transition is realized as a competition between 
two geometries (|17p . One possibility for the gravity dual to a confining theory 
at finite temperature is to take any of the metrics we described above with the 
Euclidean time periodically identified. Another option is to introduce a black 
hole horizon which completely hides the wall behind it: in Euclidean signature 
this means smoothly ending the geometry before the wall is encountered. To see 
which of those two geometries dominates the thermal ensemble one can compare 
their free energies. At low temperature the background with wall dominates 
(the hadronic phase), at large temperature the black hole (the deconfined phase). 
Out of the confining examples discussed above, an analytic solution describing 
the high temperature black hole solution is only available for Witten's black hole: 
the high temperature solution is the same as the low temperature one, just with 
the roles of spatial and time circle exchanged. Numerical black hole solutions 
have been found for the warped deformed conifold (|40|) and the N = 2* solution 
(|4ip . One pleasant feature of this transition is that the high temperature phase 
is insensitive to the details of the confinement mechanism. Many features of the 
high temperature phase are shared between all the confining theories with gravity 
dual and, in fact, also the original M = 4. Most studies of finite temperature QCD 
using the gauge-string duality use the M = 4 theory itself, as this is expected to 
be as good a model of the quark-gluon plasma as any of the confining examples. 

5 Shear and bulk viscosity 

5.1 Entropy and shear viscosity 

Special interest attaches to the ratio rj/s = 1/4-7T of shear viscosity to entropy 
density predicted by the gauge-string duality (f4~2j) . On the theoretical side, this 
calculation is appealing because of its universality (f4~3|) . and because of the con- 
jecture (|44h that rj/s > 1/4tt for any physically reasonable medium. On the ex- 
perimental side, it is interesting because elliptic flow in heavy-ion collisions seems 
to be best understood in terms of nearly ideal fluid dynamics, with < rj/s ^ 0.2, 
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a range clearly consistent with the value 1/4-7T from the gauge-string duality in 
the strong coupling limit. A recent review centered on the 77/s calculation has 
appeared in Ann. Rev. Nucl. Part. Sci. (|45p . to which we refer the reader both 
for a more thorough exposition and for a more extensive set of references. 

To understand the result rj/s = 1/4tv in a pedestrian fashion, let's start by 
considering the computation of the entropy of near-extremal D3-branes. More 
precisely, let's start with the ^cLSs-Schwarzschiid metric, ([8]), times a five-sphere 
5 5 of radius L. The horizon is at u = u^. The dt 2 part of the metric vanishes at 
this value of u, and the du 2 part becomes singular. In a more proper treatment, 
the t-u plane would be parametrized instead by Kruskal coordinates, and the 
metric would be seen to be perfectly smooth at the horizon. What matters for 
us, however, is that the dx\ + dx\ + dx 2 part of the metric is smooth at u = u^, 
and that the S 5 part of the metric is entirely independent of u. The metric on 
the horizon is just 

^horizon = L 2 u 2 L {dx\ + dx\ + dx 2 ) + L 2 dVL\ . (13) 

What this means is that the area of the horizon is just V^L^u^ Vol <S 5 , where 
Vol S 5 = 7r 3 L 5 is the volume of the five-sphere and V3 is the coordinate three- 
volume in the x\, X2, and X3 directions. Using the relation S = A^/AGio between 
the entropy S and horizon area Ah, we arrive at 

ss xr^7' <14) 

To obtain the familiar result s = |sf rec of (|16|) . one needs the relation T = u^/it 
and also the relation L 8 /Giq = 2iV 2 /7r 4 , which is specific to M = 4 super- Yang- 
Mills and AdS§ x 5 5 and originates in the formula for the tension of D3-branes. 

Next we need to understand the gravity calculation that determines 77. Suppose 
we stand far away from a near-extremal stack of D3-branes, in asymptotically flat 
space, and throw gravitons at it. When the energy ui of the gravitons is much 
less than 1/L, where L is the radius of curvature characteristic of the D3-branes, 
the absorption will be dominated by the s-wave component of the incoming flux 
of gravitons. A famous result (|2ip is that the absorption cross-section equals the 
horizon area in the extreme low-energy limit: 

^absorb = A h . (15) 

The result (|15p is intuitive: surely the s-wave cross-section would be at least 
proportional to the horizon area. But it is not inevitable that the constant 
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of proportionality should be one. This factor comes from first showing that 
the relevant linearized Einstein equations reduce to the massless Klein-Gordon 
equation in the black hole background, and then analyzing the low- frequency 
solutions of the Klein-Gordon solution with purely infalling boundary conditions 
at the horizon. 

To connect (|15p with shear viscosity, one needs three ingredients. First, we 
have Kubo's formula: 

?] = — lim — lmGu(u) , (16) 

uj^O LO 

where 

//•oo 
d 3 xj Q dt([T 12 (t,x),T 12 (0,0)]) (17) 

is the retarded Green's function of the 12 component of the stress tensor. (Any 
traceless spatial part of T mn would do as well here.) Next, we have the spectral 
representation 

ImG R (u) = ^^|(/|r 12 (0,0)K)| 2 (2^) 3 5 3 (p7-p 4 ) 

XTT[5(E f - Ei + uj)-8{E f - Ei-uj)] . (18) 

Finally, we have an application of Fermi's Golden Rule to compute the net rate 
at which gravitons are absorbed by the branes from a plane wave of gravitons 
polarized so that only h\2 is non- vanishing. The interaction of the gravitons with 
the branes is 

5 int = J d 4 xh l2 T 12 . (19) 

Far from the D3-branes, in asymptotically flat space, let's take h\2 = e _ia; ( t_a; 4) ) 
where x^ is one of the directions normal to the brane. Then the flux of gravitons 
is 

F = ^k' (20) 

The lo arises from the usual structure of fluxes for bosons: T = d </> for a 
canonically normalized complex boson cf). But h±2 isn't canonically normalized, 
due to the factor l/167rGio in front of the Einstein action in ten dimensions. The 
8-7rGio in the denominator of (|20p compensates for this, as well as the fact that 
we must take the real part of h±2 at the end of the day. Fermi's Golden Rule says 
that the net absorption rate of gravitons is 

r = v^-^\u\Ti2{^m\ 2 {^f^{p f -Pi) 

x 2vr [S(E f - Ei -co)- 5{E f - E { + io)\ 
= -2V 3 ImG R (oj). (21) 
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On the other hand, the defining relation for cr a bsorb is T = dabsorb-^"- Combining 
this relation with (|16M2ip . one arrives at 

' u^o 16ttGi 16ttGi 4vr v ' 

In the second equality of (|22p . we used the result (|15p. In the third equality, 
we used the relation S = A^/AGio. Note that we did not require the explicit 
computation (|14|) . 

The universality of the result rj/s = l/4ir may be seen as approximately parallel 
to the universality of the result cr a bsorb = ^h- More explicit checks of universality 
have been performed ([43.46.47J). Corrections from ol effects have been calculated: 
first, apparently incorrectly, in (|48|) . and later, apparently correctly, in (|49|) . Loop 
effects, corresponding to 1/N corrections in the dual gauge theory, have been 
explored as well in (|50p . and in (|5ip it was shown that the bound n/s > 1/47T 
is violated — perhaps only slightly — by a string theory construction involving D3- 
branes at an orientifold singularity The latter calculation has been checked and 
confirmed in (1521). 



5.2 Bulk viscosity 

Bulk viscosity, denoted £, quantifies energy loss due to expansion. It can be 
described in terms of the constitutive relations for relativistic hydrodynamics: 



(e + p)u m u n + pg mn - 

(23) 



pma pnb 



d a u b + d b u a - ^g a bd c u c 



Here u m = (7, jv) is the local four-velocity, and P mn = g mn + u m u n is the 
projector onto the spatial directions of the local rest frame. Indices m, n run from 
to 3, while indices a, b run from 1 to 3. Conformal invariance is summarized by 
the identity T m m = 0, and inspection of (J23J) shows that this demands the familiar 
relation p = e/3 and also £ = 0. But conformal symmetry may be broken in the 
gauge-string duality: indeed, there is a large literature studying renormalization 
group flows, and all the gravity duals of confining gauge theories are at some level 
examples of renormalization group flows. A number of authors have considered 
bulk viscosity computations in non-conformal gauge-string backgrounds ([53 1 154 1 

MJ). 

One of the main guidelines to the output of such computations is the bound 
56J) 



„* 2 U-*)- (24) 
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Computations of £ in theories where the equation of state is close to that of 
QCD have appeared in (|57|,I58[) : see also the closely related work (|59j) . where the 
equation of state is closer to that of pure Yang-Mills theory. The outcome of 
these calculations is that £/s does reach a maximum in the middle of the decon- 
finement crossover, but its maximum value is approximately bounded by 1/4-7T, 
which means that bulk viscosity is not a dominant source of entropy production 
in heavy ion collisions. There is some tension between this result and the results 
of (|6U|) , whose sum-rule calculations suggest that bulk viscosity becomes the dom- 
inant correction to inviscid hydrodynamics near freeze-out. Lattice calculations 
on pure Yang-Mills (|6ip find values of £/s as high as 0.7 very close to T c . It may 
be perilous, however, to make direct comparisons between pure Yang-Mills the- 
ory and QCD for quantities evaluated right near T c , given that one theory goes 
through a first order phase transition while the other goes through a cross-over. 

6 Expanding Plasmas 

As with perturbative methods, it is much simpler to study static finite-temperature 
configurations in the gauge-string duality than to study expanding plasmas. But 
the expansion is obviously a crucial feature of experimental conditions. A notable 
attempt to go beyond static configurations is the work (f62~j), in which an approxi- 
mate holographic dual of Bjorken flow was found. This flow has boost invariance 
in the X3 direction and translation invariance in the x\ and X2 directions. It de- 
scribes a plasma which is infinite in all three spatial directions and expands in the 
X3 direction, with a local velocity i>3 oc £3 at any fixed time t > 0. The late-time 
behavior found in ()62p is precisely what follows from inviscid hydrodynamics; 
however, a more refined analysis fj63|) shows that shear viscosity corrections enter 
with a magnitude as expected from the relation rj/s = l/4ir. This is reassuring, 
because the original derivation (|42p was for perturbations around a static configu- 
ration, contrasting with the dynamical setting of Bjorken flow. Subsequent work, 
including (64), confirms earlier suggestions, for example in (|65h . that there is a 
general correspondence between non-static black hole horizons and fluid dynam- 
ics in a dual gauge theory. The general flavor of this connection is to develop the 
dynamics of the black hole horizon in a long-wavelength expansion. The earliest 
terms in such an expansion reveal thermodynamic properties and inviscid fluid 
dynamical behavior; higher order corrections allow the extraction not only of 
the shear viscosity, but also of dynamical transport coefficients that characterize 
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the response to higher derivative deformations of uniform flow profiles. These 
transport coefficients can also be extracted from a sufficiently precise analysis of 
linearized perturbations of the AdSs-Schwarzschild solution (|66p . 

Another line of thought (|67p is to describe a finite-sized, expanding plasma in 
terms of a static black hole in global anti-de Sitter space. This proposal follows 
earlier work (|68p. and it leads to exact solutions for the expectation (T mn ) of 
the dual stress tensor. A disadvantage is that the exact solutions have perfect 
radial symmetry. They describe a spherical shell of confer mal plasma collapsing 
to a minimum radius and then re-expanding. There is no shear and no entropy 
production. However, one can perturb the exact solutions to study both hydrody- 
namical deformations and the rate of local equilibration to a hydro dynamical de- 
scription. On the dual black hole side the perturbations are quasi-normal modes 
(QNM's). As was already understood in (69), a small subset of the QNM's de- 
scribe perturbations in linearized hydrodynamics, while the others describe non- 
hydro dynamical behaviors. The rate of decay of the non-hydrodynamical modes 
is quite fast: their e-folding time is at most r e _f id ~ l/8.6T pca k, where T pea k is 
the maximum temperature of the equilibrated plasma. In (|67p it was suggested 
that, given this e-folding time, it would take a total time r ~ 0.3fm/c for the 
strongly anisotropic momentum-space distribution of the initial state to equili- 
brate to a locally thermal state under conditions similar to a gold-gold collision 
at y / s/v r /v ~ 200 GeV. A similar point was made in (63) in the context of the 
holographic dual to Bjorken flow. Related estimates have appeared in (I70ll7ip . 
The strong-coupling estimates are significantly shorter than straightforward esti- 
mates from pQCD, which — according to (|72p — are typically in the range 2.5fm/c 
or higher. 

Yet another related line of thought is to relate the formation of a black hole 
horizon in a high-energy collision in AdS§ to the production of entropy in a 
heavy-ion collision. Contributions in this area include the works (j73|) . which 
hinges on a proposed relation between gravitational shock waves in a cut-off 
version of AdS§ and a coherent pion field; (|74h . which proposes a translation of 
the whole history of a heavy- ion collision into gravitational language; (|75p . which 
provides a quantitative estimate of the total multiplicity in gold-gold collisions; 
and (|76|) . which includes a study of the metric shortly after the collision as well 
as an estimate of the thermalization time which, like the ones discussed in the 
previous paragraph, is short compared to perturbative estimates. Recently, the 
formation of a black hole in AdS§ in response to a time-dependent metric on the 
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boundary was studied numerically (|77p. The dual gauge theory equilibrates to 
an effectively thermal state in a time which is short compared to perturbative 
estimates, indicating at least qualitative agreement with the analysis of QNM's. 

7 Hard Probes and Heavy Mesons 
7.1 Hard Probes 

In QCD, the coupling runs as a function of scale. In relativistic heavy ion col- 
lisions, different sub-processes happen at different scales and hence at different 
couplings. Most of the collisions apparently lead to end-products which quickly 
thermalize and produce the expanding fireball of strongly-coupled quark-gluon 
plasma. But some of the partons collide basically head on, leading to stan- 
dard perturbative QCD scattering, which produces very energetic partons. Even 
though these energetic partons originate in the same nucleus-nucleus collision, 
they are better thought of as external probes traveling through the plasma. In 
fact these "hard probes" serve as an excellent tool to study the plasma. The 
way they interact with the plasma and the rate at which they deposit energy and 
momentum reveal detailed properties of the plasma. 

In order to get a theoretical handle on the properties of hard probes, one 
first wants to understand the idealized scenario of a probe traveling through an 
infinitely extended, spatially uniform plasma. The resulting transport properties 
can than be included in hydrodynamic codes modeling the real fireball in a heavy 
ion collision. At weak coupling, the energy loss by heavy and light quarks or 
gluons is well understood in QCD, as we partly review below. In the regime of 
the strongly coupled quark-gluon plasma, we are lacking any reliable theoretical 
tools, so getting some guidance from gauge/gravity models to the qualitative 
features of such probes is clearly desirable. 

As we discussed, in the simple example of N = 4 SYM the coupling does 
not run as a function of scale, and so one does not automatically produce hard 
partons in a high-energy collision. This point has been made most explicitly in 
the paper (|78j) . There it was shown that at zero temperature a single off-shell 
photon produces a completely round energy distribution. This is in stark contrast 
to QCD where the same process (which experimentally is realized in DIS) for an 
off-shell photon with energy far above the strong coupling scale produces a jet- 
like antenna pattern, that is energy sharply localized along the trajectories one 
would want to associate to a "back-to-back" quark/anti-quark pair produced in 
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a perturbative QCD process. In strongly coupled N = 4 SYM it simply produces 
many very soft gluons. 

In order to circumvent this issue in the study of hard probes using the gauge / 
gravity correspondence, two main routes have been pursued. One is to augment 
N = 4 SYM with additional very heavy flavors using the technique describe in 
section [3) As long as the mass is sufficiently large compared to the temperature^] 
the heavy quark will automatically serve as a hard probe of the plasma. Its energy 
loss is well described by the "trailing string" of refs ([701150]) as we will describe 
below: see also the closely related calculation of (j8ip . An earlier approach to 
energy loss from hard probes (82), loosely related to the trailing string, is based 
on the gauge-string treatment (I83p of radiation from an accelerating charge. For 
the trailing string, not just the rates, but the full energy and momentum profile 
of the wake left behind by the probe are well understood. 

Another route in the string theoretic study of hard probes is to relax the 
requirement that one must describe an entire scattering process in terms of the 
gauge-string duality. One appeals to the weakly coupled regime of QCD for 
the initial production of the hard probe and then only studies the ensuing time 
evolution of the hard probe in the strongly coupled M = 4 plasma. One example 
in this spirit is the calculation of the jet quenching parameter of (|84h . where it is 
first argued that in QCD the energy loss of a highly relativistic parton is given 
in terms of a light-like Wilson loop and then only the latter is calculated using 
the gauge gravity correspondence. Another example is the falling strings of refs 
(|85U86ll87|) where one sets up in the N = 4 SYM an initial state that corresponds 
to a sharply localized initial distribution of energy and momentum, as from an 
energetic light quark or gluon. The subsequent evolution of such a distribution 
in the M = 4 plasma can be studied completely within the framework of the 
gauge-string duality. The production of such an initial state is quite unnatural 
in strongly coupled N = 4 SYM. But such states are of course typical in weak 
coupling QCD scattering. Once produced, their subsequent evolution should be 
similar in strongly coupled QCD and strongly coupled N = 4. 

Let us describe some of the constructions in more detail. As we argued above, 

the simplest question to address is energy loss by a heavy quark. In this case 

In QCD the requirement would simply be that the quark mass be much larger than the 
temperature, but in theories with a gravity dual the meson mass is parametrically smaller than 
the quark mass due to the large binding energy. In this case the relevant criterion for whether 
a quark is heavy is whether the meson masses are much larger than the temperature. In the 
M = 4 case this for example implies that the quark mass is much larger than \f\T, not just T. 
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the problem can be entirely phrased within the N = 4 theory. Heavy flavors 
can easily be incorporated using the D7 probes of section [3l On the supergravity 
side, one is looking for a string profile that has a single endpoint on the brane 
moving along at a constant velocity v. As the quark loses energy and momentum 
at a constant rate due to the drag force experienced in the plasma, one needs 
to turn on a background electric field on the brane in order to ensure that the 
quark can maintain constant velocity; or one can work in a formal limit where 
the quark is infinitely heavy, and so maintains constant v despite a finite rate 
of energy loss. Asking for a stationary embedding of the string worldsheet of 
the form x(t,z) = vt + f(z), where x is the direction of motion the quark, z 
is the radial AdS^ coordinate and t is time, the minimal area requirement gives 
a single ODE for f{z). Furthermore, since the metric is independent of x, the 
area of the worldsheet only depends on /' and not on /, and so the equations 
of motion following from the minimal-area action for the string can simply be 
integrated. The generic solution describes a string that turns around and has 
a second endpoint on the brane/boundary: a quark/antiquark pair. There is a 
unique choice for the integration constant in /' that leads to a string with only 
one endpoint at the brane. In the IR this "trailing" or "dragging" string logarith- 
mically approaches the horizon ^| Energy and momentum flow down at a constant 
rate down the string. This corresponds to a redistribution of energy along the 
string, from parts of the string close to the endpoint (which are best thought of 
as part of the quasi-particle consisting of the quark and the surrounding gluon 
cloud) to parts of the string close to the horizon (which are best thought of as 
part of the plasma). In the case of the trailing string, the induced metric on 
the worldsheet is a 2d black hole (|89ll88|) . and so a natural choice would be to 
consider the worldsheet horizon as the boundary between the two regions, but in 
general there is some ambiguity in how to separate the string into quasi-particle 
and plasma. This is of course in complete analogy with field theory, where the 
definition of a quasi-particle as a localized lump of energy and charge suffers from 
a similar ambiguity. 

A series of papers ([901 ED EH E31 ES ESI ESI EH EH]) explains how to calculate 
the spatial distribution of energy loss starting from the trailing string. In the 
infrared, there is a good match to hydrodynamics with rj/s = 1/Att. The sonic 
boom and the diffusion wake are of comparable strength in the hydrodynamic 
regime, in line with earlier hydrodynamic calculations, in particular "Scenario 1" 



As explained in (|88[1. it crosses the horizon when described using Kruskal coordinates. 
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of ([99]) . It has been suggested (|100H99j) that a sonic boom is responsible for the 
experimental observations summarized as "jet-splitting:" namely, a broad peak 
in di-hadron histograms consistent with high-angle emission from an energetic 
parton traversing the quark-gluon plasma (jl01t lT02) . It appears, however, to be 
challenging to get a strong enough sonic boom to make this idea work (11031) . It is 
therefore interesting to note that analyses (104,105) combining the flow patterns 
calculated from the trailing string with Cooper-Frye hadronization (|106p show a 
significant jet-splitting effect from region close, but not too close, to the quark. 
Comparison of the results of (|104]I105]) with data is complicated by the facts that 
expansion of the plasma was neglected, and that the trailing string describes 
heavy quarks, which presumably make a very small contribution to untagged 
di-hadron histograms. 

The energy loss experienced by a light quark is a little more complicated. In 
QCD energy loss by an ultra-relativistic parton can be determined in the BDMPS 
framework of (110711108"]) . which describes radiative energy loss. The important 
parameter that sets the loss rate is the jet quenching parameter q = (p^_)/x, 
the average transverse momentum acquired by a parton after it has traveled a 
distance x, measured in the rest frame of the plasma. This calculation is based 
on perturbative QCD, and therefore this framework is certainly not applicable 
when discussing the energy loss of relativistic partons in strongly coupled N = 4 
SYM, whose coupling remains constant even at high energies. It was however 
argued in (84, 109) that q can be rephrased as an expectation value of a partially 
lightlike Wilson loop. This Wilson loop then can be determined in AdS/CFT 
via a minimal area calculation. Some doubts remain whether the saddle point 
identified in (|84p is the dominant contribution to the bulk path-integral, see for 
example (1110]I1H"]) . For heavy quarks, the jet quenching parameter — if defined 
as q = {p\)/x — can be directly determined in AdS/CFT ([891188]) . In fact, both 
transverse and longitudinal momentum dispersion can be determined in terms of 
small fluctuations around the trailing string described above. For a heavy quark 
with boost factor 7, the transverse momentum broadening grows as 7 1 / 2 , while 
the longitudinal momentum broadening grows with 7 5 / 2 . So formally, in the limit 
of an ultra-relativistic particle the jet-quenching parameter diverges. This is not 
in direct contradiction to the result of (I84p . as the analysis of the trailing string 
based on classical worldsheets is expected to receive quantum corrections for 7's 
that are parametrically large compared to A, whereas the calculation of q based 
on lightlike Wilson lines is best thought of being done at infinite 7 for large but 
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finite A. Last but not least, one can directly study the energy loss of a light quark 
or gluon in the strongly coupled M = 4 directly. In this case one is studying a 
falling string in the AdSs-Schwarzschild geometry. Analytic estimates of such 
falling strings from ref ()85p revealed that the stopping distance of a light quark 
is expected to scale with the energy of the quark as E 1 ^ 3 . This scaling has also 
been seen analytically and numerically in (|112l[86ll8T|) . (|57|) also revealed that 
most of the energy of the light quark is released in a burst of energy at the very 
end of the trajectory, a phenomenon that is well known as a "Bragg peak" in 
other contexts. 

7.2 Heavy Flavor and Mesons 

So far we have been discussing the properties of colored hard probes travers- 
ing the medium. For certain physical phenomena, such as J ftp suppression, it 
is also important to understand the properties of color neutral mesons moving 
through the plasma. While the high temperature phase of QCD is a theory of 
free quarks and gluons, there is some evidence from the lattice (for a review 
see (|113p ) that bound-states involving heavy flavors such as charm survive to 
temperatures above the crossover at least as peaks in the spectral density. In a 
generic finite temperature theory, there can never be sharp bound-states at finite 
temperature; thermal fluctuations always lead to a finite width. However, in the 
limit of a large number of colors that underlies the classical supergravity limit 
of gauge-string constructions, heavy flavor mesons can survive for some range of 

— 1/2 

temperatures as true bound states with a width of order N c ' that formally 
goes to zero. So the gauge-string duality is a perfect testing ground for mapping 
out properties of these heavy flavor mesons at finite temperature. 

In the study of flavors via probes branes, it has been found that a generic 
flavored gauge theory with gravitational dual exhibits a meson melting phase 
transition as a function of temperature. For temperatures below a critical value 
set by the mass of the lightest heavy flavor meson, the mesons are stable. As 
the temperature is increased they fall apart, and even heavy flavor only exists as 
free quarks. This meson melting happens in a sharp first order phase transition 
with a kink in the free energy. This phenomenon has been first observed for the 
D3/D7 system in (I114pil5|) and then has been shown to be a generic outcome of 
a probe brane analysis in (|116p . The fact that the transition is first order is an 
artifact both of the limit of a large number of colors and large 't Hooft coupling 
and it has been argued that both finite N c (|117p and finite coupling effects (|118p 
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will smooth out the transition. 

Besides the existence of this meson-melting phase transition, another inter- 
esting generic feature of mesons from probe branes is the existence of a speed 
limit. This can be seen either from the analysis of the modes living on the brane 
(I119tll20"|) . or from the analysis of moving quark/anti-quark pairs as being de- 
scribed by semi-classical strings ()121|, l79 t ri22 . 123). An important concept that 
allows one to see this effect in the dual bulk geometry is the local speed of light. 
An object moving in the spatial direction along a slice of constant radial coordi- 
nate in the black hole background given by metric (JSj) moves at the speed of light 
if its velocity is h(u) l l 2 . In other words, any object moving at a constant speed v 
along the slice is actually superluminal if it reaches down to an effective horizon 
at v 2 = h(u C h)- Using h = 1 — u^/u 4 , one obtains = n^/(l — i; 2 ) 1 / 4 . For a 
flavor brane corresponding to heavy flavors with mass m reaching down to a finite 
u m ~ ^/Xm, mesons corresponding to a small fluctuation propagating along the 
brane with a speed v experience an effective speed limit that requires u m > u^ ff jf| 
A similar speed limit is in effect for large spin mesons or a quark-antiquark dipole 
that is held up by an external field. In that case, the meson or dipole is described 
by a semi-classical string in the bulk that reaches down to some minimum turn- 
ing point uq which is set by the separation between quark and antiquark. Again, 
the speed limit requires uq > u^ ff . While the details depend on the setup, the 
fact that the speed limit compares the typical energy scale of the excitation not 
to T but to T/(l — u 2 ) 1 / 4 is universal. Similar scaling can be observed in the 
fluctuation analyses of (|89U88p . 

8 Concluding remarks 

The dream of describing strong interactions in terms of string theory is an old 
one — older than QCD itself. The gauge-string duality, as the modern incarna- 
tion of that dream, still falls short in important respects: asymptotic freedom is 
poorly understood, and gauge-string constructions typically lead to a parametric 
mismatch between the mass gap and the flux tube tension. These defects are 
due to the almost universal reliance on the classical supergravity approximation, 
corresponding to the large N, large gy M N limit of the dual gauge theory. And 
yet, as we have reviewed, constructions using classical gravity and probe branes 

3 A phenomenological application has been proposed in (p-24): a flavorless meson, like the 
J/V'j may produce an observable peak in the photon spectrum if its limiting in-medium velocity 
is subluminal. 



Pedestrian's Guide 31 

provide an excellent venue for calculations ranging from meson spectra to trans- 
port coefficients. Finite temperature is elegantly included in these calculations 
in terms of a black hole horizon in the fifth dimension. Making the comparisons 
between such calculations and QCD systematic is a persistent difficulty. Despite 
this difficulty, gauge-string treatments of confinement and of finite-temperature 
non-abelian plasmas are quite valuable because they complement insights from 
more standard quantum field theoretic treatments. For example, the descriptions 
of confinement are elegant and geometrical; the connection to hydrodynamics is 
relatively simple; and the interplay between hard probes and a thermal medium 
is rich and explicit. 

Attempts to forecast the future of a field of theoretical physics usually are 
doomed to failure. However, we look forward to future developments in the con- 
nections between AdS/CFT and strong nuclear interactions in several directions, 
including: studies of thermalization in terms of dual black holes; some better uni- 
fication of the descriptions of energy loss; progress, if only incremental, toward 
finite 't Hooft coupling; and better integration of the "top-down" constructions 
of confining duals with "bottom-up" models of holographic QCD. 
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